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Abstract 



The transform considered in the paper integrates a function sup- 
ported in the unit disk on the plane over all circles centered at the 
boundary of this disk. Such circular Radon transform arises in several 
contemporary imaging techniques, as well as in other applications. As 
it is common for transforms of Radon type, its range has infinite co- 



■ dimension in standard function spaces. Range descriptions for such 
transforms are known to be very important for computed tomogra- 

iy-j \ phy, for instance when dealing with incomplete data, error correction, 

■ and other issues. A complete range description for the circular Radon 
. transform is obtained. Range conditions include the recently found 

set of moment type conditions, which happens to be incomplete, as 
well as the rest of conditions that have less standard form. In order 
to explain the procedure better, a similar (non-standard) treatment of 
the range conditions is described first for the usual Radon transform 
^ ' on the plane. 

1 Introduction 



The following "circular" Radon transform, which is the main object of study 
in this article, arises in several applications, including the newly developing 
thermoacoustic tomography and its sibling optoacoustic tomography (e.g., 
EH |HU E3 EH1 EH ED] ) , as well as radar, sonar, and other applications 
[JOl HH1 EH E2j • It has also been considered in relation to some problems of 
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approximation theory, mathematical physics, and other areas |T| [21 EH HS1 

[2H1 EH EH1 EH] 

Let f(x) be a continuous function on M d , d > 2. 
Definition 1. JTie circular Radon transform of f is defined as 

Rf(p,p)= [ f(y)da{y), 

where da(y) is the surface area on the sphere \y — p\ = p centered at p G M d . 

In this definition we do not restrict the set of centers p or radii r. It is 
clear, however, that this mapping is overdetermined, since the dimension of 
pairs (p,r) is d+ 1, while the function / depends on d variables only. This 
(as well as the tomographic motivation) suggests to restrict the set of centers 
to a set (hypersurface) S C M d , while not imposing any restrictions on the 
radii. This restricted transform is denoted by Rs'- 

Rsffap) = Rf(p,p)\ P es- 

In this paper we will be dealing with the planar case only, i.e. the di- 
mension d will be equal to 2. Due to tomographic applications, where S is 
the set of locations of transducers [2H1 ETJ EH UH] , we will be from now on 
looking at the specific case when S is the unit circle |x| = 1 in the plane. 

There are many questions one can ask concerning the circular transform 
Rs'- its injectivity, inversion formulas, stability of inversion, range descrip- 
tion, etc. Experience of computerized tomography shows (e.g., jUJEH]) that 
all these questions are of importance. Although none of them has been re- 
solved completely for Rs, significant developments have occurred recently 

(e.g., [H [21 EE El H21 [HI Eini Em ESI ESI EH EU E21 Eni EH [HHl [EH]) The goal 

of this article is to describe the range of Rs in the two-dimensional case, 
with S being the unit circle. Moreover, we will be dealing with functions 
supported inside the circle S only. The properties of the operator Rs (e.g., 
stability of the inversion, its FIO properties, etc.) deteriorate on functions 
with supports extending outside S (e.g., [21 EH EUJ and remarks in the last 
section of this article). However, in tomographic applications one normally 
deals with functions supported inside S only [23 EH EH IZDJ - 

As it has already been mentioned, the range of Rs has infinite co-dimension 
(e.g., in spaces of smooth functions, see details below) and thus infinitely 
many range conditions appear. It seems to be a rather standard situation for 
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various types of Radon transforms that range conditions split into two types, 
one of which is usually easier to discover, while another "half" is harder to 
come by. For instance, it took about a decade to find the complete range 
description for the so called exponential Radon transform arising in SPECT 
(single photon emission computed tomography) El HUH EH EH] • For a more 
general attenuated transform arising in SPECT, it took twice as much time 
to move from a partial set of range conditions |Hl to the complete set 
1531 . In the circular partial set of such conditions was discovered 

recently [HS] • It happens to be incomplete, and the goal of this text is to find 
the complete one. 

One might ask why is it important to know the range conditions. Such 
conditions have been used extensively in tomography (as well as in radiation 
therapy planning, e.g. [HI EH [30] ) for various purposes: completing incomplete 
data, detecting and correcting measurement errors and hardware imperfec- 
tions, recovering unknown attenuation, etc. jMIITll^lOIIlliniloTllloTllMl 
155] . Thus, as soon as a new Radon type transform arises in an application, 
a quest for the range description begins. 

In order to explain our approach, we start in the next section with treating 
a toy example of the standard Radon transform on the plane, where the 
range conditions are well known (e.g., [JJl E3 CHI EH 123 EH ESI)- Our 
approach, however, is different from the standard ones and naturally leads 
to the considerations of the circular transform in the rest of the paper. 

2 The case of the planar Radon transform 

In this section we will approach in a somewhat non-standard way the issue 
of the range description for the standard Radon transform on the plane. 
Consider a compactly supported smooth function f(x) on the plane and its 
Radon transform 



where s G f , w 6 S 1 is a unit vector in M 2 , and dl is the arc length measure 
on the line x • u — s. We want to describe the range of this transform, say 
on the space C^°(]R 2 ). Such a description is well known (e.g., [TBI IT71 ITS} IT§| 
[23 EH ES] , or any other book or survey on Radon transforms or computed 
tomography): 




(1) 
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Theorem 2. A function g belongs to the range of the Radon transform on 
Cq° if and only if the following conditions are satisfied: 

1. ge C^iS 1 x R), 

oo 

2. for any k G Z + the k-th moment Gk{oj) = j s k g(u, s)ds is the restric- 

— oo 

tion to the unit circle S 1 of a homogeneous polynomial of ui of degree 
k, 

3. g(u,s) = g{-u),-s). 

We would like to look at this result from a little bit different prospective, 
which will allow us to do a similar thing in the case of the circular Radon 
transform. 

In order to do so, let us expand g(u, s) into the Fourier series with respect 
to the polar angle ip (i.e., uj = (cos if), sin ■?/>)) 

oo 

g{u,s)= 9n(s)e in t (2) 

n=— oo 

We can now reformulate the last theorem in the following a little bit 

strange way: 

Theorem 3. A function g belongs to the range of the Radon transform on 
if and only if the following conditions are satisfied: 

1. ge C^iS 1 x R), 

oo 

2. for any n, the Mellin transform Mg n (o) = j s a ~ 1 g n (s)ds of the n- 

o 

th Fourier coefficient g n of g vanishes at any pole a of the function 

r(^fW), 

3. g(u,s) = g(-uj,-s). 

Since the only difference in the statements of these two theorems is in the 
conditions 2, let us check that these conditions mean the same thing in both 
cases. Indeed, let us expand g(w,s) into Fourier series (J2J) with respect to 
if). Representing e m ^ as the homogeneous polynomial {u)\ + i(sign n)^)'"' of 
uj of degree \n\, and noticing that uj\ + uj\ = 1 on the unit circle, one easily 
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concludes that the condition 2 in Theorem|2]is equivalent to the following: the 
k-th moment J s k g n (s)ds of the n-th Fourier coefficient vanishes for integers 

M 

< k < \n\ such that k — n is even. 

Let us now look at the condition 2 in Theorem El still using the same 
Fourier expansion. Notice that when k — \n\ is a negative even integer, 
Mg n (a) is one-half of the moment of order k = o — 1 of g n . Taking into 
account that T( a+1 ~^ ) = r( fc+2 ~^ ) has poles exactly when k — \n\ is a 
negative even integer, we see that conditions 2 in both theorems are in fact 
saying the same thing. 

One can now ask the question, why should one disguise in the statement of 
Theorem El negative integers as poles of Gamma- function and usual moments 
as values of Mellin transforms? The answer is that in the less invariant 
and thus more complex situation of the circular Radon transform, one can 
formulate a range description in the spirit of Theorem El albeit it is unclear 
how to get an analog of the version given in Theorem 

As a warm-up, let is derive the condition 2 in Theorem El directly, without 
relying on the version given in the preceding theorem. This is in fact an easy 
by-product of the A. Cormack's inversion procedure, see e.g. jH3 Section 
II. 2]. Indeed, if we write down the original function f(x) in polar coordinates 
r(cos0,sin0) and expand into the Fourier series with respect to the polar 
angle 

oo 

/(r(cos0,sin0))= £ f n {r)e m t (3) 

n=— oo 

then the Fourier coefficients /„ and g n of the original and of its Radon trans- 
form are related as follows |U1 formula (2.17) and further]: 

"(r/„(r)) W = (4) 



B r , 



where 



T(s)2~ s 

B n (s) = const—- : — — : — ,, . , (5) 

r((s + 1 + H)/2)r((s + 1 - \n\)/2) v ; 

Thus, condition 2 of Theorem El guarantees that the function M(rf n (r))(s) 
does not develop singularities (which it cannot do for a Q^-function /) at 
zeros of B n (s). It is not that hard now to prove also sufficiency in the 
theorem, applying Cormack's inversion procedure to g satisfying conditions 
1-3. However, we are not going to do so, since in the next sections we will 
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devote ourselves to doing similar thing in the more complicated situation of 
the circular Radon transform. 



3 The circular Radon transform. Formula- 
tion of the main result 

Let us recall the notion of Hankel transform (e.g., JTj). For a function h(r) 
on M + , one defines its Hankel transform of an integer order n as follows: 

oo 

(H n h)(a) = J J n (ar)h(r)rdr, (6) 
o 

where the standard notation J n is used for Bessel functions of the first kind. 

Let, as in the Introduction, R$ be the circular Radon transform on the 
plane that integrates functions compactly supported inside the unit disk 
D over all circles \x — p\ = p with centers p located on the unit circle 
S = {p\\p\ = 1)- Since this transform commutes with rotations about 
the origin, the Fourier series expansion with respect to the polar angle par- 
tially diagonalizes the operator, and thus the n-th Fourier coefficient g n (p) of 
g = Rsf will depend on the n-th coefficient f n of the original / only. It was 
shown in jHT] that the following relation between these coefficients holds: 

g n (p) = 2npH {J n H n {f n }}. (7) 

For the reader's convenience, we will provide the brief derivation from |5T] . 
Considering a single harmonic / = f n {r)& in ^ and using polar coordinates, 
one obtains 



2tt 

intf> 



gn(p) = J rf n (r)dr J 5 [(r 2 + 1 — 2rcos0) 1 ^ 2 — p\ e 
o o 



Thus, the computation boils down to evaluating the integral 



2tt 



/ = J 8 [(r 2 + 1 - 2rcos0) 1/2 - p] e 
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Using the standard identity 

S(pf - p) = p I J {p'z)J (pz)zdz 



and the identity that is easy to obtain from one of the addition formulas, e.g. 
from formula (4.10.6)] 

2-Kj n {az)J n {bz) = J J [z(a 2 + b 2 - 2a6cos0) 1/2 ] e - in 
o 

one arrives from (JSJ) to (JZJ). 

Since Hankel transforms are involutive, it is easy to invert (JJJ) and get 
Norton's inversion formulae joTJ 

/« = ^Wn{ j- }• (9) 

Now one can clearly see analogies with the case of the Radon transform, 
where zeros of Bessel functions should probably introduce some range con- 
ditions. This happens to be correct and leads to the main result of this 
article: 

Theorem 4. In order for the function g(p,p) on S 1 x R to be representable 
as Rsf with f e Cq 3 (D), it is necessary and sufficient that the following 
conditions are satisfied: 

1. geCf^x (0,2)). 

oo 

2. For any n, the 2k-th moment J p 2k g n (p)dp of the n-th Fourier coeffi- 

o 

cient of g vanishes for integers < k < \n\. (Equivalently, the 2k-th 

oo 

moment J p 2k g(p, p)dp is the restriction to the unit circle S of a (non- 
o 

homogeneous) polynomial of p of degree at most k.) 

oo 

3. For any n G Z, function TCo{g n (p) / p}(cr) = J Jo(o'p)g n (p)dp vanishes 

o 

at any zero a ^ of Bessel function J n . (Equivalently, the nth Fourier 
coefficient with respect to p 6 S 1 of the "Bessel moment" G a (p) = 

oo 

J J (ap)g(p, p)dp vanishes if a ^ is a zero of Bessel function J n .) 
o 
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4 Proof of the main result 



Let us start with proving necessity, which is rather straightforward. Indeed, 
the necessity of condition 1 is obvious. Let us prove the second condition. In 
fact, it has already been established in [HE]- Let us repeat for completeness 
its simple proof. Let k be an integer. Consider the moment of order 2k of g: 



p 2k g{p,p)dp= \x -p\ 2k f{x)dx = {\x\ 2 - 2x ■ p + l) k f(x)dx (10) 



R2 R2 

(we have taken into account that \p\ = 1). We see that the resulting expres- 
sion is the restriction to S 1 of a (non-homogeneous) polynomial of degree k 
in variable p. Expanding into Fourier series with respect to the polar angle 
of p, we see that the nth harmonic g n contributes the following homogeneous 
polynomial of degree \n\ in the variable p: 



P 2k g n (p)dp 



Here as before p = (cosip, smtp). Thus, for |n| > k, this term must vanish, 
which gives necessity of condition 2. We will return to a discussion of this 
condition below to add a new twist to it. 

Necessity of condition 3 follows immediately from Norton's formula (JJ2J), 
which implies in particular that 

n {g n (p)/p} = 2nJ n H n {f n }. 

Since both functions J n and TC n {f n } are entire, Tlo{g n (p) / p} vanishes when- 
ever J n does. 

Remark 5. The reader might ask why in the third condition of the Theorem 
we do not take into account the zero root of J n , which in fact has order n, 
while non-zero roots are all simple. The reason is that the condition 2 already 
guarantees that a = is zero of order 2n of TCo{g n (p) / p} (twice higher than 
that of J n ). Indeed, due to evenness of J , function T~£o{g n (p)/ p}{o~) is also 
even. Thus, all odd order derivatives at a = vanish. The known Taylor 
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expansion of Jo at zero leads to the formula 



DC 



. ^ ( — ]) m /r r\ 2m f 

n {g n (p)/p}(a) = Y: (g) / r 2m g n (r)dr. 
m 

We see now that the moment condition 2 guarantees that a = is zero of 
order 2n of H {g n / p}(a). 

Let us move to the harder part, proving sufficiency. Assume a function g 
satisfies conditions of the theorem and is supported in S x (e, 2 — e) for some 
positive e. We will show that then g = R$f for some / G C™(D e ), where D t 
is the disk \x\ < e in the plane. 

Due to Norton's formulas, it is natural to expect the proof to go along 
the following lines: expand g into the Fourier series g = Y g m {p)& 1 with 

m 

respect to the angle variable ip, then use to construct a function / and 
then show that / is of an appropriate function class and that its circular 
Radon transform is equal to g. This is what we are going to do, with a 
small caveat that instead of constructing / itself, we will construct its two- 
dimensional Fourier transform. Besides, we will start considering the partial 
sums of the series h n — Yl gm(p)e im< ^- But first, we need to get some simple 

\m\<n 

estimates from below for the Bessel function of the first kind J n . 

Lemma 6. On the entire complex plane except for a disk So centered at the 
origin and a countable number of disks Sk of radii 71/6 centered at points 
7f(k + one has 

n \Imz\ 

1^(^)1 > OO (11) 



Proof: Let us split the complex plane into three parts by a circle So 
of a radius R (to be chosen later) centered at the origin and a planar strip 
{z = x + iy\ \y\ < a}, as follows: part I consists of points z satisfying \z\ > R 
and \Imz\ > a; part II consists of points such that \z\ > R and \Imz\ < a; 
part III is the interior of So, i.e. \z\ < R. It is clearly sufficient to prove 
the estimate (jll)) in the first two parts: outside and inside the strip. Using 
the parity property of J n , it suffices to consider only the right half plane 
Rez > 0. 
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The Bessel function of the first kind J n (z) has the following known asymp- 
totic representation in the sector | argz\ < it — 5 (e.g., formula (4.8.5)] or 
[HI formula (5.11.6)]): 

J„(: ) = \flf z cos(z - f - f)(l + 0{\z\-*)) 

(12) 



sin{z-?f-l)(^ + 0{\z\ 



Let us start estimating in the first part of the complex plane, i.e. where 
\Im z\ > a and \z\ > R for sufficiently large a and R (and, as we have agreed, 
Rez > 0). There, due to boundedness of tan z in this region, one concludes 



that = cos z (0(\z\ *)), and thus (fT2*|) implies 

z 

J n {z) = \Pfcos(z - ^ - + oor 1 )), 

V ixz l 4 

which in turn for sufficiently large a, R leads to 



n p \Imz\ 

\J n {z)\ > (13) 



In the second part of the plane (right half of the strip), due to boundedness 
of sin z we have 



2 



Uz) = \ - 




ITT) 7T 

cos(z - — --)(! + 0(\z\- 2 )) + OQzt 1 ) 



TTZ I 2 4' 

Consider the system of non-intersecting circles Sk with centers at Zk = \ + 
kir + ^Y + j and radii equal to |. Then outside these circles | cos(2; — z y 1 — 1)| > 
C and 

^^^^(i+od^r 1 )). 

\z\ 



This implies that for a suitably chosen and sufficiently large R, inside of the 
strip and outside the circles Sk, we have 

n p \Imz\ 

\J n (z)\ > ^j=^ (14) 



for | z | > .R. This proves the statement of the lemma. □ 
Let us now return to our task: consider the function g and the partial 
sums h„ of its Fourier series. 
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Lemma 7. 1. If g((f),p) = J2 m 9m(p)e im ^ satisfies conditions of Theo- 
rem^ and is supported in S x (e, 2 — e), then each partial sum h n = 

T J \m\<n9m{py imi ' does SO. 

2. For any n, h n is representable as Rsf n for a function f n G C^(D e ). 

Proof of the lemma. The first statement of the lemma is obvious. 

Due to i), it is sufficient to prove the second statement for a single 
term g = g n (p)e tn ^ '■ As it was just mentioned, we will reconstruct the 
Fourier transform F of the function /. In order to do this, we will use 
the standard relation between Fourier and Hankel transforms. Let as before 
f(x) = f n (r)e m< ^, where r = \x\ and are polar coordinates on R 2 . Then the 
Fourier transform F(£) of / at points of the form £ = au, where a G C and 
u = (cos if), sinip) G M 2 can be written up to a constant factor as follows: 

F(au)=H n (f n )(a)e i ^. (15) 

(e.g., [TTJ end of Section 14.1]). If we knew that g = Rsf, then according to 
(0) this would mean that 

F(au) = F(a)e m ^ = J_ ^o(gn(p)/p)(a) ^ 

27T Jn(V) 

Let us now take this formula (|16|) as the definition of F(au). Due to the 
standard parity property of Bessel functions, such F is a correctly defined 
function of auj for o ^ (i.e., F(au) = F((—a)(—u))). We would like to 
show that it is the Fourier transform of a function / G C^°(D t ). Let us prove 
first that F belongs to the Schwartz space iS(M 2 ). In order to do so, we need 
to show its smoothness with respect to the angular variable if), smoothness 
and fast decay with all derivatives in the radial variable a, as well as that no 
singularity arises at the origin, which in principle could, due to usage of polar 
coordinates. Smoothness with respect to the angular variable is obvious, due 
to (|16|) . Let us deal with the more complex issue of smoothness and decay 
with respect to a. First of all, taking into account that g n (p) is supported 
inside (0, 2), and due to the standard 2D Paley- Wiener theorem, we conclude 
that u(a) = Ti (g n (p)/p) is an entire function that satisfies for any iV the 
estimate 

\u(a)\ < C N (1 + H)-V 2 - £)|/mff| . (17) 

According to the range conditions 2 and 3 of the Theorem, this function 
vanishes at all zeros of Bessel function J n (a) at least to the order of the 
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corresponding zero. This means, that function F(a) = — - is entire. 

2nJ n {a) 

Let us show that it belongs to a Paley- Wiener class. 

Indeed, 7i(g n (p)/p) is an entire function with Paley- Wiener estimate (JT7J). 
Due to the estimate from below for J n (jllj) given in Lemma El we conclude 
that F(ouj) is an entire function of Paley- Wiener class in the radial directions, 
uniformly with respect to the polar angle. Namely, 

\F(a)\ < CW(1 + \a\)- N e^ Ima \. (18) 

Indeed, outside the family of circles Sk the estimate (fTT|) together with 
(fT7j) give the Paley- Wiener estimate (fT£J) we need. Inside these circles, appli- 
cation of the maximum principle finishes the job. Smoothness with respect 
to the polar angle is obvious. Thus, the only thing one needs to establish 
to verify that F belongs to the Schwartz class is that F is smooth at the 
origin. This, however, is the standard question in the Radon transform the- 
ory, the answer to which is well known (e.g., jT7| p. 108-109], [111113, 123 
Ch. 1, proof of Theorem 2.4]). Namely, one needs to establish that for any 
non- negative integer k, the kth radial (i.e., with respect to cr) derivative of 
F(au) at the origin is a homogeneous polynomial of order k with respect to 
uj. So, let us check that this condition is satisfied in our situation. First of 
all, the parity of the function F is the same as of n. Thus, we do not need 
to worry about the derivatives F^ with k — n odd, since they are zero 
automatically. Due to the special single-harmonic form of F, we only need 

(k) 

to check that F& |o-=o = for k < \n\ with k — n even. This, however, as we 
have discussed already in Remark follows from the moment conditions 2 
of the Theorem. 

Due to the smoothness that we have just established and Paley- Wiener 
estimates, F € 5(1R 2 ). Thus, F = f for some / G iS(M 2 ). It remains to 
show that / is supported inside the disk D e . Consider the usual Radon 
transform Hf(s,(f>) of /. According to the standard Fourier-slice theorem 
[EH H3 UHl EH1 123 133! , the one-dimensional Fourier transform (denoted by a 
"hat" ) from the variable s to a gives (up to a fixed constant factor) the values 
TZf(a, if)) = F(au)), if as before u = (cosip, smif>). Here 71, as before, denotes 
the standard Radon transform in the plane. Since functions F(au) of cr, as 
we have just discussed, are uniformly with respect to to of a Paley- Wiener 
class, this implies that lZf(s,u) has uniformly with respect to to bounded 
support in \s\ < 1 — e. Now the "hole theorem" (231331 (which is applicable 
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to functions of the Schwartz class), implies that / is supported in D e . 

The last step is to show that Rgf = g = g n {r)e tn ^. This, however, 
immediately follows from comparing formulas (jl6j) and (|7J), which finishes 
the proof of the main Lemma □ 

Let us now return to the proof of Theorem 0] We have proven so far that 
any partial sum h n of the Fourier series for g belongs to the range of the 
operator Rs acting on smooth functions supported inside the disk D e . The 
function g itself is the limit of h n in C^(S x (e, 2 — e)). The only thing that 
remains to be proven is that the range is closed in an appropriate topology. 
Microlocal analysis can help with this. 

Consider Rs as an operator acting from functions defined on the open 
unit disk D to functions defined on the open cylinder Q = S x (0,2). As 
such, it is a Fourier integral operator j^H 123 EH] ■ If Rs is the dual operator, 
then E = R s Rs is an elliptic pseudo-differential operator of order —1 [2*T| 
Theorem l] 1 . 

Lemma 8. The continuous linear operator E : Hq(D e ) i— > Hf oc (D) has zero 
kernel and closed image. 

Proof of the lemma. Since E = R s Rs, the kernel of this operator 
coincides with the kernel of Rs acting on Hl(D e ). Since 5* is closed, it is 
known that Rs has no compactly supported functions in its kernel |2] 
(this also follows from analytic ellipticity of E and Theorem 8.5.6 of [27], see 
also Lemma 4.4 in [2 ). Thus, the statement about the kernel is proven and 
we only need to prove the closedness of the range. 

Let P be a properly supported pseudo-differential parametrix of order 1 
for E [HH]. Then PE = I + B, where B is an infinitely smoothing operator on 
D. Consider the operator II that acts as the composition of restriction to D e 
and then orthogonal projection onto Hq(D^) in H 2 (D e ). On Hq(D e ) one has 
UPE = I + K, where K is a compact operator on Hq{D^). Notice that the 
operator IIP is continuous from the Frechet space Hf oc (D) to Hq(D c ). Due 
to the Fredholm structure of the operator UPE = I + K acting on Hq(D e ), 
its kernel is finite-dimensional. Let M C Hq(D e ) be a closed subspace of 
finite codimension complementary to the kernel, so / + K is injective on M 
and has closed range. Then one can find a bounded operator A in Hq^D^) 
such that A(I + K) acts as identity on M. Thus, the operator AUP provides 

1 Bolkcr's injective immersion condition [23 that is needed for validity of this result, is 
satisfied here, as shown in the proof of Lemma 4.3 in 0. 
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a continuous left inverse to £ : M h Hf oc (D). This shows that the range of 
E on M is closed in Hf oc (D). On the other hand, the total range of E differs 
only by a finite dimension from the one on M. Thus, it is also closed. □ 
We can now finish the proof of the theorem. Indeed, the last lemma shows 
that the function R s g, being in the closure of the range, is in fact in the range, 
and thus can be represented as Ef with some / G Hq(D € ). In other words, 
R s (Rsf — g) = 0. Since the kernel of R s on compactly supported functions 
is orthogonal to the range of Rs, we conclude that Rsf — g = 0. Since 
Ef = R s g is smooth, due to ellipticity of E we conclude that / is smooth as 
well. This concludes the proof of the theorem. □ 

5 Remarks and acknowledgments 

We would like to finish with some remarks. 

• It should be possible to prove that the operator Rs in the situation 
considered in the text is semi-Fredholm between appropriate Sobolev 
spaces (analogously to the properties of the standard and attenuated 
Radon transforms, e.g. [2U IS])- This would eliminate the necessity 
of the closedness of the range discussion in the end of the proof of 
Theorem 0] 

Such a statement could probably be proven either by using FIO tech- 
niques, or by controlling dependence on n of the constant C and of the 
radius of the circle So in Lemma El The former approach would be 
better, being more general. 

• Proving compactness of support of function / in Lemma [TJ we used the 
standard Radon transform and the "hole theorem." Instead, one could 
probably use the fact that Fourier transform of / is, by construction, 
a Paley- Wiener class CR-function on the three-dimensional variety of 
points au in C 2 and then use an appropriate mandatory analytic ex- 
tension theorem in the spirit of [o^j . 

• We considered the situation most natural for tomographic imaging, 
when the functions to reconstruct are supported inside the aperture 
curve S. What happens when the supports of functions extend outside 
the circle SI It is known that compactly supported [2] (or even belong- 
ing to L p with sufficiently small p [1]) functions can still be uniquely 
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reconstructed. Necessity of the range conditions we derived apparently 
still holds and they are still sufficient for finite Fourier series. However, 
many things do go wrong in this case. Our proof of the closedness of 
the range fails (in particular, since the Bolker's condition for the corre- 
sponding FIO does not hold anymore, which was also the main hurdle 
in proving the results of j2J). Moreover, the range will not be closed 
anymore. Indeed, reconstruction will become unstable, since due to 
standard microlocal reasons [321 HHH HOI EH ED] , some parts of the wave 
front set of the function outside S will not be stably recoverable. This 
means, in particular, that non-smooth functions can have smooth cir- 
cular Radon images. This, in turn implies that the range is not closed 
in the spaces under consideration, and so sufficiency of the range con- 
ditions should fail. We are not sure what kind of range description, if 
any, could work in this situation. By the way, the nice backprojection 
type inversion formulas available in odd dimensions JE] also fail for 
such functions. 

• It would be interesting to understand range conditions in the case of a 
closed curve S different from a circle. Since our method uses rotational 
invariance, it is not directly applicable to this situation. 

• Our result is stated and proven in 2D only. It is possible that a similar 
approach might work in higher dimensions. As we have been notified by 
D. Finch, he and Rakesh have recently obtained by different methods 
some range descriptions in 3D [To] . 

• The result of the paper was presented at the Fully Three-Dimensional 
Image Reconstruction Meeting in Radiology and Nuclear Medicine July 
6-9, 2005 in Salt Lake City, Utah. 

This work was supported in part by the NSF Grants DMS 9971674 and 
0002195. The authors thank the NSF for this support. Any opinions, find- 
ings, and conclusions or recommendations expressed in this paper are those 
of the authors and do not necessarily reflect the views of the National Science 
Foundation. 
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